MULTISPECIES VIRIAL EXPANSIONS 

SABINE JANSEN, STEPHEN J. TATE, DIMITRIOS TSAGKAROGIANNIS, AND DANIEL UELTSCHI 

Abstract. We study the virial expansion of mixtures of countably many different types 
of particles. The main tool is the Lagrange-Good inversion formula, which has other ap- 
plications such as counting coloured trees or studying probability generating functions in 
multi-type branching processes. We prove that the virial expansion converges absolutely 
in a domain of small densities. In addition, we establish that the virial coefficients can 
^ be expressed in terms of two-connected graphs. 
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1. Introduction 



(~| One of the central questions of statistical mechanics is the calculation of the thermody- 

C^ namic properties of a given fluid starting from the forces between the molecules. A major 

(— I contribution in this direction is the work of Mayer and his collaborators |May37| IMA371 

" ' IHM38I IMH38J . This was further developed by Born-Fuchs |BF38j and Uhlenbeck-Kahn 

|UK38j ■ All of these papers concern the theory of non-ideal gases for which the pressure is 

written as a power series in the activity or density. The convergence of such power series was 

investigated much later, during the 1960's, by Groeneveld, Lebowitz, Penrose and Ruelle. 

Parallel to the study of monoatomic gases, multi-component systems were also investigated, 

K^ for example in |May39| for the case of a two-component system and later in jFuc42| for a 

mixture with an arbitrary (but finite) number of different components. Although briefly 

mentioned already in |Fuc42| , the complete study of the convergence of the activity expan- 

^N^ sion comes later in |BL64) for mixtures of finitely many components. 

■ ' The generalisation of virial expansions and approximate equations of state from the 

f~^ monoatomic gas to mixtures is not straightforward. The van der Waals equation for binary 

(^ mixtures of hard spheres, for example, comes in different versions, distinguished by different 

^-H "mixing rules" for obtaining effective parameters of the mixture |HL70) : see also |LR64j 

^ for the Percus-Yevick and virial equations of state. Going from binary mixtures to multi- 

component systems is not easy either - in Fuchs' words, "for the most intricate parts of the 
calculation even the theory of the two-component system does not give any indication as to 
jrt the results for the many-component system (e.g. the reduction of the cluster integrals to 

the irreducible cluster integrals)" |Fuc42| . Even less trivial is the extension from finitely to 
infinitely many types of particles, in particular for estimating the domain of convergence. 

The present work deals with the virial expansion for countably many types of particles, 
addressing in particular the problem of convergence. Our motivation is two-fold. First, 
mixtures as studied in physical chemistry are of interest on their own. Second, multi- 
component systems arise as effective models as in "the consideration of certain phenomena 
connected with the order-disorder transition in alloys" |Fuc42] . or in the treatment of a 
monoatomic gas as a mixture of "droplets" (groups of particles close in space) |IIil56[ Chapter 
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5.27]. A droplet can comprise arbitrarily many particles, and there can be arbitrarily many 
droplet sizes - this is why we are interested in mixtures with infinitely many species. 

In the first part of this article, we are given a pressure function that depends on countably 
many fugacity parameters. This yields countably many densities (that are also functions of 
fugacities) and the goal is to write the pressure as function of densities. In the case of one 
parameter, this can be done using Lagrange inversion |LP64j . We note that this analytical 
method is also a standard tool in combinatorics iBLL98[ Chapter 3] . In the case of multi- 
component systems, the key ingredient is Good's generalisation of the Lagrange inversion 
to several variables IGoo d60J . The Lagrange-Good inversion has attracted attention in a 
variety of contexts [AbdOSl fBru83l [Ges87l IEM94| . Faris has recently noticed its relevance for 
the virial expansion |Farl2j . Our main focus is on the convergence of the expansion. Despite 
recent renewed interest on the single-species virial expansion |Janl21 IPT12[ ITatelSl IMP13| , 
this question does not seem to have been addressed before in the case of infinitely many 
species. Our main result is the existence of a non-trivial domain of convergence for the 
expansion of the pressure as function of densities. A novel feature is to use Lagrange-Good 
for proving convergence. 

The second part of this article deals with a gas of classical particles with two-body 
interactions. There are many species of particles. Under some assumptions, the virial 
coefficients are given in terms of two-connected graphs (irreducible cluster integrals). Early 
derivations for systems with one, two, or finitely many components can be found in |BF38[ 
|May39|rFuc42) . Thanks to the work of Mayer a systematic connection with combinatorics has 
been initiated, linking the enumeration of connected and two-connected graphs with cluster 
and virial expansions. This connection has been further developed in the work of Leroux and 
collaborators |Ler04j . leading to modern proofs for the expression of virial coefficients. This 
was generalised by Faris to the case of many species |Farl2| . The present article has some 
overlap, but there are several differences. In particular, we formulate sufficient conditions 
on the interactions that guarantee the convergence of the virial expansion. 

We would like to emphasise that our results have relevance beyond statistical mechan- 
ics. Indeed, our main result can be formulated as an inverse function theorem for functions 



between Banach spaces that are not necessarily Frechet-differentiable (Section 2.2). In ad- 
dition, our result can be applied in the original context of the Lagrange-Good inversion 
formula: Good motivated his work by stochastic branching processes and combinatorics of 
coloured trees |Good60| IGood65J . Recursive properties of trees lead to functional equa- 
tions between generating functions. When combined with the inversion formula they yield 
expressions of probabilities or tree cardinalities as contour integrals; Good explicitly com- 
puted some of those integrals. Our result yields bounds for cardinalities without having 
to compute the integrals, which may come in handy when explicit computations prove too 
complicated. 

The article is organised as follows. In Section |2] we give the general setting of the virial 
expansion in the context of formal power series. The main theorem proposes sufficient 
conditions under which the virial power series is absolutely convergent. This is done by 
deriving explicit bounds on the virial coefficients via the Lagrange-Good inversion formula. 
In Section [3] we prove a dissymmetry theorem for coloured weighted graphs and deduce that 
the virial coefficients are given by the two-connected graphs; this result holds whenever the 
weights satisfy a block factorisation property. Finally, in Section HI we consider a mixture 
of rigid molecules and, using the results on the convergence of the cluster expansion given 
in [Uel041 IPUOQj , we show that the mixture meets the conditions of Section [2] 
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2. General virial expansions 



2.1. Setting &: results. Let z = {zi,Z2, ■ ■ ■) denote a sequence of complex numbers. Con- 
sider the formal series 

p{z)=J2Kn)z'', (2.1) 

n 

where the sum is over all multi-indices n — (ni, n2, . . . ), n^ G N, with finitely-many non-zero 
entries and X^ieN ^i — 1- We assume that for all i G N, the coefficient &(0, ...,1,0,...) of z^ 
in p{z) is non-zero; this is the only condition needed for Lagrange inversion techniques. In 
Sections [3] and l4] we will need the additional assumption that the coefficients are normalized 
as 6(0, . . . , 1, 0, . . .) = 1, as is the case in most apphcations in Statistical Mechanics. 

We use the notation z'" = J|,- z"* . In statistical physics, Zi represents the "activity" (or 
"fugacity") of the species i, and p{z) is the pressure of the system with many species. A 
physically relevant quantity is the density pi of the species i, whose definition is 

Pz[Z) ^ Z, — (Z). (2.2) 

We do not suppose yet that the series for p{z) is convergent, and the equation above should 
be understood in the sense of formal series. To be precise, Pi{z) is the formal series with 
coefficients nih{n). 

In statistical physics, the virial expansion is the expansion of the pressure in powers of 
the density. Accordingly, we define the coefficients c{n) by the equation 

P{z)=^c{n)p{z)'^. (2.3) 

n 

Here also, we use the notation p{z)'^ = Y\iPi{z)^K Let us check that the coefficients c{n) 
are well-defined in the sense of formal series. Observe that [z^]p{z)'^ ^ only if n < fc, 
i.e., if rii < ki for all i; the notation [z^]p{z)'^ refers to the coefficient of z^ in the formal 
power series p{z)'^. Then {p{z)'^)n is a "summa ble family" of formal series in the sense of 
Def. 2.2 in |EM94) . and the coefficients of z'^ in K^ satisfy 



b{k) = [z''] J2 <n)p{zr - Y. c(n)[^^]p(^)", (2.4) 



i<fc 



the latter sum being finite. Then Eq. (2.4 1 can be inverted recursively, and c(fci, . . . , /c^, 0, . . . ) 
can be expressed in terms of b{n) and c(ni, . . . , n^, 0, . . . ) with Ui < ki, i — 1, ...,£— 1, and 
nn < ki. 

Our goal is to control the convergence of the virial expansion, assuming convergence of the 
series p{z). In the following the statement "| log/(2;)| < a" means that there is a A{z) G C 
such that f{z) = exp{A{z)) and 1^(2) | < a, i.e., the precise choice of the branch of the 
logarithm is irrelevant. 

Theorem 2.1. Assume that there exist < r^ < i?,; and a^ > 0, i = 1, 2, . . . , such that 

• p{z) converges absolutely in the polydisk D — {z (£ C | Vi € N : \zi\ < Ri}. 
'1 ^P f \ 



• 



< ai for all i > 1 and all z £ D. 



dzi 

* / rr— < 00 and > „ < 00. 



r,a^ 



Ri 
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Then there exists a constant C < oo (which depends on the r.i, Ri, at, but not on n) such 
that 

\c{n)\ <C sup \p{z)\Y[( — r 
ZED ^^^ r, I 



(2.5) 



i>l 



The estimate for cin) guarantees convergence of the series ^^ c(ri)p^ for ah p in a 
polydisk 



'-{ 



D'= peC" IVieN: \p,\ <ne- 



iSN 



e ' 



< oo 



}■ 



(2.6) 



Theorem 1 2.1 1 is proved in Section [273} 

We can also address the foUowing question: Consider the functions Zi{p) obtained by 



inverting (2.2); for given p e D' , does z{p) belong to D, so that p{z{p)) is given by an 



absolutely convergent series? The following result provides a partial answer, as it guarantees 
convergence when p belongs to a smaller domain. 



Theorem 2.2. Under the same assumptions as in Theorem \2.1\ we have 

,z{p)'^ _ -n-/ e"' \ni+fc. 



[P" 



<C 



n(^)' 



The constant C is the same as in Theorem |2.1[ and the proof is similar. It can be found 
at the end of Section [23] Theorem [2.2[ implies that for all p <E D' , 



\Z^{P)\ < C\p, 



J>1 



\pj\ 



<CR,\{[1- 



\P3\ 



(2.7) 



so that z{p) G D for p small enough. The inequality (2.7) is the analogue of the bound 



|pfc(^)| < l^fcl exp(afc), valid for z € D under the assumptions of Theorem 2.1 



2.2. The point of view of the inverse function theorem. The formulation of The- 
orem 2.1 is geared towards the virial expansion in statistical mechanics. The theorem in 



itself is, however, purely analytic. In this section we rephrase it as a type of inverse function 
theorem and discuss its relation to traditional inverse function theorems. 

Let (^Fk{w)Y be a family of power series in the complex variables Wj [j e N) such that 
Ffe(O) 7^ 0, for all fc e N. The goal is to invert the system of equations u^. = WkFk{w). On 
the level of formal power series, the inversion is always possible: there is a unique family of 
power series Gk{u), fc S N, such that the inverse is given by Wk{u) = UkGk{u), i.e., for all 
A: € N, we have Wk{u)Fk{w{u)) = Uk, as an identity of formal power series. 

Theorem 2.3. Assume that there exist < ri < Ri and ai > 0, i = 1,2, ... , such that 

• The series Fk{w), k E N, converge absolutely in the polydisk D — {w £ C^ | Vz € 
N: \w,\<R,}. 

• I \ogFi{w)\ < ai for all i > 1 and all w E D. 

1 / — - < oo and > — -— < oo . 

V Ri ''-^ Ri 

i>l * * i>l * 

Then there exists a constant C < oo (which depends on the ri, Ri, ai, but not on n) such 
that for all fc e N, 

.^ / pO-i \ rii 

|K]GfcH|<Ce'^'= n( . 

Z>1 
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The proof is similar to the proofs of Theorems |2.1| and |2.2[ In order to better understand 
the analytic structure of the theorem, it is convenient to introduce Banach spaces of complex 
sequences. To simplify matters, suppose that the first two assumptions of the theorem hold 
with afc = ak, where a > is some fc-independent constant, and Rk = exp(— afc). Situations 
of this type occur in the context of cluster expansions. Choose r^ := k~^ exp(— afc). Define 
the weighted ^-'^-norms ||i(;||)3^g :— "^j^^-^ k^ e^p{qak)\wk\- Let us define 

Bp,q{£) ■■= {w eC^^ : \\w\\p^g < e}. (2.8) 

As a consequence of Theorem |2.3| we get the well-defined functions 

/ : Bo,i(l) ^ Bo,o(l) g ■■ Biaie) ^ Bi^i{Cc,e) C So,i(l) 

W ^ {wkFk{w)) ^^^^ U^ {ukGk{u))^^^ (2.9) 

Then for any u G i?4^2(£), we have 



f{a{u)) =u. 



(2.10) 



Here e e (0, 1), c^ is such that for all x € (0, e), (1 — x)~^ < c^ exp(j:), and we assume that 
e is small enough so that Cc^ exp(e)e < 1. 

It is instructive to compare this with traditional inverse function theorems: suppose that 
/, considered as map from the Banach space with norm || • ||4,i to the space with norm 
II • ||4_2, was Frechet-differentiable with invertible derivative in a neighbourhood of the origin. 
Then we would get the existence of open neighbourhoods of the origin such that / is a 
bijection between these neighbourhoods |Zei95[ Chapter 4.10]. Our result yields a weaker 
conclusion: we have a bijection between i34_2(e) and g(i?4_2(e)), but in general the latter set 



needs not be open with respect to || • ||4,i. The reason is that Theorem 2.3 operates under 
conditions that are weaker than those of the inverse function theorem: for infinitely many 
variables, the existence of continuous partial derivatives df/dzk does not imply that / is 
Frechet-differentiable — we do not know whether the Jacobi matrix {dfi/dz}.) represents a 
bounded operator. Our condition exp(— afc) < \fk{z)lzk\ < exp(afc) replaces the traditional 
condition that the derivative and its inverse are bounded operators between Banach spaces. 

2.3. Lagrange-Good inversion & bounds of virial coefficients. The Lagrange-Good 
inversion formula gives explicit expressions for c{n) and [p'^]z^ , which can then be estimated. 
We use Eq. (4.5) of |Ges87| to get 

[p"]$(^(p)) = [z-Mz) j^ det M{z). (2.11) 

\dz) 

Here, we used the notation (g^) — HilaF) '• Let N be the largest index i such that 
n^ ^ 0. M{z) is the iV X A^ matrix 

Af(^)= % + ^^ . (2.12) 

^ a^ ^ i<jj<7v 

We employ the formula above with $(2:) — p{z) and $(2) = z^ . In |Ges87| . the for- 
mula has been proved for finitely many species; a proof for infinitely many species is given 
in |EM94j (Theorem 4). Note, however, that we can apply the finitely many species version 
in our context because we only need it for n with finitely many non-zero entries. 

In order to estimate the coefficients of the right side, we use Cauchy's formula and get 



upper bounds on the various terms. We start with the determinant in (2.11) 
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Lemma 2.4. Under the assumptions of Theorem \2.1[ there exists a constant C < oo such 
that for all z with \zi\ — ri, and all N , we have 



\<leiM{z)\ <C. 



Proof. We have the expansion 



deiM[z)= ^ det 



JC{l,...,A'} 



dp 



(2.13) 



Let Ui be non-zero numbers to be determined later. We use det M — det DMD ^ with D 
the diagonal matrix with entriesu"^^ in the diagonal, and we get 



detM(2;)= ^ det 

JC{1,...,N} 



Uj d T dp 



Ui ozj ozi 



(2.14) 



3 / ijej 



From Hadamard's inequality, we get the upper bound 



dzi dzj 



|2 ^ 
,2 



1/2 



(2.15) 



|detM(^)|< Y. nlEl^^l' 
Jc{i,...M}ieJ^jeJ 

By Cauchy's formula, choosing the contour around z^ with radius Ri — r^, and using the 

dp 
dzj ■ 

j_ /.log^(i«,u;) 
27ri 



bound on the logarithm of ^, we get 



d , 9p 
-—log 75— (2 
ozi ozj 



dzj 
(w - ZiY 



-dw 



< 



Ri - ri 



(2.16) 



>W ^. 



where z^ Ms the vector z without the Zi term. Then 



detM(z)|< Y. n 

JC{l,...,A'}ieJ 



Ui[Ri 



T^(^"^ 



^a^ 
"j 



jeJ 



1/2- 



N 

< exp 



1 + 



u,{R,-r,)\j^^'^^''^ 



(E^^^ 



E;7715^(E^?« 



UiiRi — ri) \^-^ ^ ^ 



1/2- 



1/2- 



(2.17) 



Choosing Uj = ^rj/Rj, the expression above is finite by the assumptions of the lemma. D 



Proof of Theorem \2.1\ We use the Lagrange-Good inversion (2.111 and Cauchy's formula 
and we write 






'^^\my 



det M{z), 



(2.18) 



with the contours being circles of radii r^ around th e ori gin. It follows from the assumptions 



that Itj^I > e "' for all i and all z £ D. Theorem 



2.1 



then follows from Lemma 



2.4 



D 
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Proof of Theorem \E^ We use Eq. ( 2.11 ) with $(2) = z'' / p{z)'' . We then repeat the proof 
of Theorem 2.1 without the term p{z). This yields the bound 



W 



Apf 



\W 



ni-\-ki 



■ dot M{z) 



1>1 ' 



(2.19) 



n 



3. Connected and two-connected graphs 

We now make a further assumption on the series p{z) introduced in ( |2.1[ ): it is given by 
the (weighted) exponential generating function of coloured graphs. This choice is motivated 
by applications to statistical mechanics, which are discussed in Section |4] If the weights 
satisfy a certain block factorisation property, the virial coefficients c{n) can be expressed 
using two-connected graphs. We therefore introduce these coloured graphs in some detail. 

Recall that a graph g is a pair g = {V{g),£{g)) with V{g) a set and £{g) C {{i,j} C 
V{g) : i y^ jY The elements of V{g) are called "vertices" and the elements of £{g) "edges". 
A graph g' is a subgraph of 5 iiV{g') C V{g) and £{g') C S{g), in which case we write g' C g. 
A graph g is connected if |V((?)| > 1 and for every partition V{g) = Vi U V2, there exists 
i e Vi and j G V2 such that {«, j} G £{g)- A graph g is two-connected if \V{g)\ > 2 and for 
every i € V(5), the subgraph g \ {i}, obtained by removing i and all edges containing i, is 
connected. 

A coloured graph is a pair (g, k) where g is a graph and and fc is a function V{g) -^ N 
that assigns the colour fc^ g N to each i £ V{g). Coloured connected graphs are pairs {g,k) 
with g a connected graph, coloured two-connected graphs are pairs (<?, fc) with g a two- 
connected graph. The results of this section can be formulated in the general framework 
of labelled coloured combinatorial species f |MN93| . |EM94[ Sections 2 and 3]), but for the 
reader's convenience we present the results in a self-contained way. 

Let w{g, fc) be a weight function on coloured graphs. We assume that it is invariant 
under relabellings that preserve the colour: let cr : V — ^ V be a bijection with the property 
that fco.(j) = ki for all i € V, and ga the graph with vertices V{ga) = V{g) and edges 
£{g„) = {{cr{i),a{j)} : {i,j} £ £{g)}, then w{ga,k) = w{g,k). The weighted exponential 
generating function for connected graphs is defined by 



^(^'-i:^ i: 



Zki ■■■Zk„ ^ ■w{g,k) 

n>l '" fc=(/ci,...,fc„) geCr, 

z"- 



(3.1) 



geC[n\ 



In the first line we denoted C„ the set of connected graphs with vertex set {1, . . . , n}. In the 
second line the sum is over multi- indices n with finitely many entries, and C[n\ denotes the 
set of coloured graphs g = {5,fc'"'') with vertices 1,2, . . . , |n| and fc'"' the colouring such 
that the first ni vertices have colour 1, the vertices ni -f 1, . . . ,7ii + 712 have colour 2, etc. 
We shall refer to fc'"' as the canonical colouring. The second expression for C{z) is more 
elegant but the first expression turns out to be more practical. These formula should still 
be understood as formal series. 
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We suppose that w{g, k) factorises with respect to the block decomposition of g. RecaU 
that an articulation point of 5 is a vertex i e V such that the subgraph g\{i} is disconnected. 
(A two-connected graph is a connected graph without articulation point.) A block is a two- 
connected subgraph g' of g that is maximal, i.e., if g" is any two-connected subgraph such 
that g' C g" C g, then g" = g' . Let {gi, . . . ,gm} be the block decomposition of g, i.e., the 
set of blocks of g. It has the following properties, (i) The blocks induce a partition of the 
edge sets: £{g) = U^^f (g^) with £(5^)0 S{gj) = when i ^ j. (ii) Each articulation point 
belongs to more than one V{gi), other vertices belong to exactly one V{gi). 





(. _ _ ^^1 


_ ' " "" ^ \ 


;' fc::;^^ J\ 


'.•\ Block 1 / / 


1 / Block 3\/ j 


\ \. / ' 
^ \ / / 

^ - Xi " "' _-— — 


'^^Q^^r^^^^' 


"^1 cVrr; 


7\"''' 


' \Biocir? 


//^ Articulation 


V N^^ y 
s ^*^^ 


/ points 




/ 



Figure 1. The block decomposition of a connected graph. 



Theorem 3.1. Assume that p{z) ~ C{z) as above, that the weight function w{g, k) satisfies 
w{g, k) — I when g has size n ^ 1, and 



'{9,k) ^Ww{9i,k\y, .) 



i=l 

when g has size n>2, where {gi, . . . ,gm} is the block decomposition of g and k\ , , is the 
restriction of the colouring k : V{g) -^ N to V{gi). Then 

fc>l n:\n\>2 geBln] 

where B[n] consists of the two-connected coloured graphs g — {g, fc'"') with canonical colour- 
ing fct"! and vertex set V{g) = {1, 2, . . . , \n\}. 

The proof is given at the end of this section. It uses the dissymmetry theorem of combi- 
natorial structures, following |BLL98[ rLer04j . This is rather straightforward but we need to 
add colours to all objects. 

If fc g N denotes a colour, a fc-rooted graph is a triplet {g,k,i) where g is a graph with 
finite vertex set V{g), k e N^^^\ i G V{g), with the property that the root i has colour 
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ki = k. The weighted exponential generating function of fc-rooted connected graphs is 



1 

~\ / J ^f^l ■ ■ ■ ^^71 I l^ ■ "'^ "'J I / J 

n>l ' ■ fc=(fci,. ..,/£„) g&Cr, 



C"'iz)^Y^^ V Zk,...Zk\{i:h = k}\y^w{g,k) 



^V^^-fc— r V" w{g) (3.2) 

•■^ — ' n! -"^ — ' 

n geC[n] 

dC 

Let C, C'^ be the sets of connected (resp. fc-rooted connected) coloured graphs with vertex 
set of the form {1, . . . ,n}, n € N, and set and C := VJk&iC'^ . The associated exponential 
generating function is C*{z) = X]fcGN^*'^(-^)- ^® define B, B* and their exponential gener- 
ating functions B{z), B'{z) in a similar way, replacing "connected" by "two-connected". 

Next, we describe the composition of connected and two-connected graphs. The set 
B{{C'''}k&i) consists of coloured two-connected graphs whose vertices contain a rooted 
connected graph with the appropriate colour for the root. More precisely, an element g € 
B[{C*''}k&i) of size n is a triple 7 = (fc, 7, (7i)iGV(7)) consisting of 

• A colour assignment k G N". 

• A two-connected graph 7 with vertex set V(7) C {1, . . . ,n}, |V(7)| > 2. 

• A family {ji)iev{-y) of connected graphs 7^ such that i e V(7i), and the vertex sets 
form a partition {1, ... ,71} = Uigv(7)"^'(7i)- Note that (7^, fc|v('yi),i) is a fc,;-rooted 
coloured connected graph. 

With each 7 we associate the connected graph g — 5(7) with vertices 1, . . . , n and edge set 
f (7) U (Uigv(7)'^(7i))- We assign to the composite structure 7 the weight ^(5(7), fc) of the 
underlying connected colored graph. We also introduce B* {{€*''} ken) , which is as above but 
with the additional choice of a root in V(7). 



Lemma 3.2. Under the assumptions of Theoreni \3.1\ the weighted exponential generating 
functions of B{{C*^}k£n) and B' {{C^} ken) satisfy 

V -.Zk, . . . ZkMg. k) = B{{C"'iz)}keN)- 

^ — ' n! 

•7eB({C"=}(ceN) 

Here, n = n{g) is the size of g, k = k{g) is the colour assignment of the vertices {1, . . . , n} 
and g — .9(7) is the induced connected graph on the set of vertices {1, . . . ,n}. Further, we 
also have 



^-Zk, . . . ZkM9,k) ^ B'{{C"'iz)}kef 
n\ 

TGB'({C"=}fceN) 






The proofs are tedious but immediate: one sums over all components of g, and uses the 
multinomial theorem so that the elements of the partition become independent. This is 
possible because of the factorisation property of the weights, 

^i'(ff(7),fc) =«^(7,fc|v(7)) X n "'(7*)'s|v(7,))- (3-3) 

iGV(7) 

For a proof in the context of labelled coloured combinatorial species (but without weights) , 
see |MN931 Proposition 1.4] and |EM941 Propositionl.3]. 

Next, we state the dissymmetry theorem for coloured graphs. Here, A + B denotes the 
disjoint union where the elements of A and B are distinct by definition. 
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Theorem 3.3. We have 

C + B'{{C"'}keN)^C'+B{{C"'}keN) 

in the sense that there is a size and weight preserving bijection between C + B*{{C*''}kef>]) 
andC'+BiiC^jkeN)- 

Proof. Let 7 h- > (7(7) be the mapping assigning a connected graph 5(7) to the compos- 
ite structure 7 as described above. The map extends naturally to two maps (j) : C + 
B'{{C'''}keN) -> C and ip '■ C* + l^{{C*''}kGn) -^ C to coloured connected graphs. It is 
sufficient to prove that for every g = {g, k) S C, the preimages 4>~^{{g}) and ^p^^{{g}) have 
the same cardinality. Indeed, if this is true, there is a bijection between the preimages, and 
the collection of bijections defines a bijection between C+;B*({C*'^}fcgN) and C*+B {{€*''} keN) 
with the required properties. 

Let g = {g,k) e C. If g has size one, then 0~^({g}) ~ {g} C C and t/j'^dg}) = 
{(g, k, 1)} C C , and both preimages have cardinality 1. If g has size n > 2, let {gi, . . . , gm} 
be the block decomposition of g. A moment of reflection shows that the preimage oi g under 
(j) consists of 

• The graph itself, g G C. 

• The composite structures (fc,7, (7i)iev(7):'') ^ B'{{C'''}keN) where 7 is one of the 
blocks gi,...,grn oi g, r e V(7), and (7^) are uniquely determined by g and the 
choice of 7. 

It follows that the preimage under (j) has cardinality 

m 
l + Y.\^(9^)\■ (3.4) 

The preimage under ip has cardinality 

n + TO. (3.5) 

The first "n" corresponds to the number of ways to choose the root when the preimage is 
in C*, and to is the number of composite structures 7 G B* {{€'''} ken) with g{j) — g. 
There remains to show that the block decomposition of any g G C„ satisfies 

m 

l + ^\V{gt)\=n + m, (3.6) 

or equivalently 

J2{\V{g,)\-l)=n-l. (3.7) 

i=l 

This can be seen by induction. This clearly holds when m — \ and V((?i) = {1, . . . ,n} 
(this corresponds to g being two-connected). Now suppose that 5 G C„ has to blocks of size 
ni , . . . , n„i . Consider the bipartite graph t whose vertex set consists of the blocks (71 , . . . , g„^ 
and articulation points ai,...,af of g, and edges {{(7^,0^} : aj G gi}. The graph t is 
known pLL98. Section 4.3] to be a tree and called the block cut-point tree oi g. Let v G V{t) 
be a leaf of t. Then d is a vertex belonging to exactly one edge and must be a block v = gi 
containing exactly one articulation point a of g. 

Thus there is a block containing precisely one articulation point of g, and without loss 
of generality we take this to be the TOth block. We remove from g all edges of the block g^ 
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and all vertices of g™, except the articulation point a. We now have a graph with m — 1 
blocks and so we have X^ta ('^i — 1) = n — n„j by induction. Therefore, 



/(jii — 1) = n — n„j + n„j — I ~ n — 1. (3.8) 






Proof of Theorem \3.1\ Lemma [3. 2 1 and the dissymmetry theorem imply that the exponential 
generating functions satisfy 

C{z) + B'{{C'Hz)}kefi) = C'{z) + i?({C"=(^)}feeN). 0.9) 

We have p{z) = C{z) and pk{z) — C'^{z), and the theorem follows. D 

We conclude this section with two remarks. The first remark is that under the assump- 



tions of Theorem 3.1 we also have a formula for the expansion of the chemical potential 



log Zfc in terms of the density, 



log Zk^ log Pk{z) - -^{p{z)), B{p)= XI ^ XI ^(^)- 

^'^ n:\n\>2 ' geB[n\ 



This follows from the relation 



Pk{z) = C'^iz) = 2;fcexp(-— (p(2;; 



(3.10) 



(3.11) 



|FarlO[ Section 3.2]; dkB is the generating function for two-connected graphs whose root is 
a "ghost" of colour k. 

The second remark is that Theorem 13.11 is not limited to connected and biconnected 
graphs, but holds for pairs of combinatorial structures with a similar composition structure 
- as Leroux puts it, for "various tree- like structures" |Ler04j . A well-known example is the 
dissymmetry theorem for trees |BLL98[ Chapter 4] , which for coloured trees reads as follows: 
Let t{z) be the exponential generating function for coloured labelled trees (unrooted), and 
E2{p) ■= \ Z]fe Pk+UkKt PkPi = hiJ2k PkY the exponential generating function of coloured 
sets of cardinality 2 (no weights). We have dkE2{p) = J2k Pk^ ^^'^ Tk{z) := Zk{dkt){z) is the 
generating function for coloured trees with a root of colour k. The analogue of Eq. ( |3.11 ) 
reads 

Tk{z) = Zfcexp(X^j(^)) (3.12) 

j 
(see for example [FarlOi Section 3.1]) and the analogue of the virial expansion is 



k<l 

2 (3-13) 



t{z) =Y^Tk{z) -Y^'-^ -Y^Tk{z)T,{ 

k k 

fc k 

Adding {J^k'^f^i^))'^ t° both sides, we obtain 

tiz) + (y^Tkiz)) =Y,Tk{z)+E2{{Tk{z))keN) (3.i4) 



which, noting that X]fe(^fe-^2)(p) — (X^fe Pfc)^i is nothing else but the dissymmetry theorem 
for (coloured) trees, compare |BLL98[ Chapter 4]. 
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In the previous example, our results add little to what is already known for trees. It 
would be interesting to know, however, which information can be obtained for coloured 
trees with colour-dependent weights and constraints - for example, the modified choice 
E2{p) '■= '^k<i PkPt corresponds to trees in which adjacent vertices have always different 
colours. This is particularly interesting because Good's original motivation for his multi- 
variable version of the Lagrange inversion came from branching processes in probability and 
combinatorics of trees |Good60[ IGood65j . 

4. Classical gas of rigid molecules 

We now describe a physical system that fits the theory of Sections [2] and |3] It consists of 
a gas of molecules that are assumed to be rigid. Let A C K'' be the domain, which we take 
as a cube in M'* with periodic boundary conditions. We let V denote its volume. A molecule 
is represented by 

X = {k,x,(j)), (4.1) 

where k & N denotes the species, x & A denotes the position, and (/i e $ := S"^^^ ={</)€ 
M.'^ I II 011 = 1} denotes the orientation. Interactions are given by a function U{Xi,X2) that 
takes values in M U {-|-oo}. Let 

C(Xl,X2)=e-^(^-^^) -I. (4.2) 

We take periodic boundary conditions, i.e., we assume that if A = [0, L]'', and y — y' E 
LZ'^, then U{{x,k,(l)),{y,£,ip)) = U{{x,k,(j)),{y' ,£,ip))- We make three assumptions on the 
interactions. The first one is about symmetries, the second one is the stability condition 
that ensures the existence of the thermodynamic limit, and the last one implies that we 
consider a regime of low density or high temperatures. 

Assumption 1. The potential function U satisfies 

• Symmetry: U{Xi,X2) = U{X2,Xi). 

• Translation invariance: If X + a denotes the molecule translated by a d A, i.e., with 
position X + a, then U{Xi + a, X2 + a) = U{Xi, X2). 

• Rotation invariance: If RX denotes the molecule rotated by the orthogonal matrix 
R, i.e., with orientation R(j), then U{RXi, RX2) ~ U{Xi, X2). 

The partition function of the system is 
^a(-2) = ^ -f ^ Zfc^ ...Zfc, / da;i...dx„/ d0i . . . d(^„ expj - ^ U{Xi,Xj)\. 

n>0 fceN" "^^^ "^^ l<i<j<7i 

(4.3) 

The term n = is understood to be equal to 1, and d(j) is the unique rotationally invariant 
measure on $ with J^ d(f> = 1. Given a graph g e C„, we define the weight function WA{g, k) 
to be 

w\{g,k)^—l dxi...dxn di^i . . . d0„ TT ({Xi,Xj), (4.4) 

U-j}e£{g) 

where 

C(X„^,)=e-^(^-^^)-l. (4.5) 

The empty product is set to be equal to I, so that graphs of size I have weight V^^ /^ dx J^ d(f) 
1. By a standard cluster expansion, or by the exponential formula of combinatorial struc- 
tures, we have 



Za{z) =expl V^— Y^ Zk,...Zk„ Y^ ■WA{g,k] 



I .^ , ^^ -. ..„ ^^ ..w, , , (4-6) 

^ n>\ fcGN" geC„ 
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— pVpA 



We then 



(4.7) 



The partition function is related to the finite-volume pressure by Za 
define 

Pa(z) = C{z), 

where C{z) is the exponential generating function of connected graphs given in (3.1) with 
the weights WA{g, k) in ( 4.4 ). The goal is to show that the assumptions of Theorem 2.1 hold 
true uniformly in the volume V . 

Lemma 4.1. Under As sumption\R the weight function of {4-4^ satisfies the block factori- 
sation 



WA 



(5,fc) = n'^A(5i,fc)- 



It is not too hard to check that factorisation holds when the graph is cut at any articulation 
point, and the lemma follows. Next, the stability condition. 



Assumption 2. There exists a nonnegative constant b such that for all n and all Xi, 
we have 

n 
l<i<j<n i—1 

In addition, we also assume that for all X and Y of species k and i, we have 

|l + C(^,>')| < c^™"('=''^). 



,Xn, 



(4.8) 



(4.9) 



The next and last assumption is the "Kotecky-Preiss criterion" that guarantees that the 
interactions and the weights are small. 

Assumption 3. There exist positive numbers Ri, R2, . . . and a constant a such that for all 
X = {k,x,(f)), 

Y^ Efc- e('^+3^)'=' f dx' f d0'|C(X,X')| < ak. 

where b is given in Assumption[^ In addition, we also assume that 

J2 Rk' e(''+3'')'^-' < 00. 



(4.10) 



(4.11) 



fc'GN 



Theorem 4.2. Let pa{z) = C{z) and suppose that Assumptions\m^hold true. Then 
(a) pa{z) converges absolutely in the polydisc D ~ {z E C^ : \zi\ < i?i Vi G N}. 
'"•^ < ak for all z ^ D and all k E N. 



(b) 



log ^(2;) 



The main consequence of this theorem is that Theorem 2.1 applies, hence the existence 
of a domain of densities with absolute convergence of the virial expansion. 

Proof. The setting of [PU09| applies by taking X = N x A x $ for the space of "polymers" 
and by defining the measure on X as 



I ■d^ = y^Zk I dx I i 



(4.12) 



The conditions of Theorem 2.1 in f PUOQj are fulfilled and it implies that for any {ki,xi, (pi), 
we have 

J2n Y. \zkJ...\z,jY\'"iS:k)\<{e'"'' ~l)e""'\ (4.13) 

n>2 fc2,...,fe„eN g<£C„ 

for all z e D. 
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For (b) we need to control the logarithm of the derivative of pA- It is not entirely 
straightforward as we need both lower and upper bounds for g^PA- We have 



dph 19 1 1 9Za 

dzk V dzk V Za{z) dzk 



(4.14) 



From the definition (4.3 1 of the partition function, we get 



dzk 



(2) = y V 7 TTT V Zfei . . . Zfc„ / dx2 . . . dx,„ / 



fc2,...,fc„>l 



n 



(4.15) 



i=2 



2<i<j<n 



We set X = (fc,0,0). The formula holds because of translation and rotation invariance, 
and because /d(/>i = 1. We observe that ^|^ is a partition function where each molecule 
Xj gets the extra factor Q-^'^^'^i) . We can again perform a cluster expansion or use the 



exponential formula of combinatorial structures. It is indeed convergent thanks to (4.9 1. We 
get 



dzk 



(2;) = yexp-^ V" — V" Zki...Zk„ da;i...da;„/ 



n 



-u{x,x 



sec„ {»j}e£(g) '' 



(4.16) 



This allows us to combine it with the cluster expansion of Za in (4.14) and we get 



dpA 
dzk 



(2;) = exp<^ V" — V" 2;fcj...Zfe„ / da;i...dx„ / 
141 ""■ fci,.tCeN -^A" ^ 



(n(l + C(X,X,))-l) ^ n C(^-^^) 
Next we use the identity 

n n—1 

[](i + C(^,x,))-i= [[](i + c(x,Xj))-i](i + c(x,x„))+c(x,x„). 

It allows to prove by induction that 



(4.17) 



(4.18) 






(4.19) 
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The integrand of (4.171 is then less than 



dx2 . . . dx„ I d(l>2 ■ ■ . A4>n\ ^ W C{Xi,Xj) 

9eC„ {ij}G£(s) 



A"-i J*"-i 



(4.20) 



< ^ |zfeje('^+3'')'=i / da;i / d^i|C(X,Xi)| 

,. ^,x, ^A J$ 



< afc. 

We bounded the bracket by e(°+2^)'=i , see Eq. (2.9) in JPU09) . The last inequality follows 
from Assumption [3] D 
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